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ABSTRACT 

Recent high redshift surveys for 21-cm absorption in damped Lyman-a absorption systems (DLAs) 
take the number of pubhshed searches at Zabs > 2 to 25, the same number as at Zabs < 2, ahhough the 
detection rate at high redshift remains significantly lower (20% cf. 60%). Using the known properties 
of the DLAs to estimate the unknown profile widths of the 21-cm non-detections and including the 
limits via a survival analysis, we show that the mean spin temperature/covering factor degeneracy at 
high redshift is, on average, double that of the low redshift sample. This value is significantly lower 
than the previous factor of eight for the spin temperatures and is about the same factor as in the 
angular diameter distance ratios between the low and high redshift samples. That is, without the 
need for the several pivotal assumptions, which lead to an evolution in the spin temperature, we show 
that the observed distribution of 21-cm detections in DLAs can be accounted for by the geometry 
effects of an expanding Universe. That is, as yet there is no evidence of the spin temperature of gas 
rich galaxies evolving with redshift. 

Subject headings: quasars: absorption lines — radio lines: galaxies — galaxies: star formation — 
galaxies: evolution — galaxies: ISM — cosmology: observations 



1. INTRODUCTION 

The 21-cm transition of hydrogen (Hi) traces the cool 
component of neutral gas throughout the Universe, the 
raw material fuelling star formation. In emission, the low 
probability of the 21-cm spin-flip transition compounded 
by the inverse square law, renders this essentially unde- 
tectable at z > 0.2. In absorption, however, this sen- 
sitivity limit can be avoided since the line strength is 
independent of distance, depending only upon the col- 
umn density of the absorbing gas and the flux from the 
background emission source. Nevertheless, detections re- 
main rare with only 77 Hi 21-cm absorption systems 
known at z > 0.1 — 35 associated with the quasar host 
galaxy (zabs = 2qso)^j in addition to 42 systems lo- 
cated at some intermediate redshift along the sight-line 
to the quasar (zabs < -zqso)^- These intervening ab- 
sorbers occur in known or candidate^ damped Lyman-a 
absorption systems (DLAs, defined by their large neu- 
tral hydrogen column densities, A^hi > 2 x 10^° cm~^). 
Currently t here arc 1500 DLAs known (jCurran et al.l 
120021: iNoterdaeme et al.ll2009t) , which could contain over 
80% of the neutral g as mass density of the Universe 
(jProchaska et al.ll2005i ). 

The velocity integrated optical depth of the 21-cm 

absorption, t = — In^l — is related to the total 

neutral hydrogen column density via (jWolfe fc Burbidg^ 
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^ Summarised inlCurran fc Whiting (2010), with the addition of 
three new associated absorbers reported in Curran ct al. ( 2011a j). 

^ Summarised in Curran (2010), w ith the addition of those re- 
ported bylCurran ct al. (2011b); Sria nand et al.l (120121 ). 

^ Usually in Mgll absorption systems at z^bs ~ 1-7, thus not 
being sufficiently redshifted to be detected in the optical band by 
ground-based telescopes. 



where Tgpin [K] is the mean harmonic spin temperature 
of the gas, a is the depth of the line (or r.m.s. noise in 
the case of a non-detection) and 5 & / the flux density 
and covering factor of the background continuum source, 
respectively. This latter parameter accounts for the pos- 
sibility that the absorbing gas may not fully cover all of 
the background radio emission'* and so all of the mea- 
sured flux, 5, may not be intercepted, decreasing the 
sensitivity of the observation. 

In both the associated and intervening cases the de- 
tections occur overwhelmingly at z < 1. For the asso- 
ciated systerns, iVn T is generally not known, although 
ICurran et al.l (120081 ) attribute the lack of absorption at 
high redshift to the traditional optical selection of tar- 
gets biasing searches towards the most luminous objects 
at these large distances, where the ultra-violet luminosity 
ionises the neutral gas. 

In the optically thin regime (ct ^ f-S),^ Equ. [T] 
reduces to Nm 1.823 x 10*^^^ I^dv, thus giv- 
ing a degeneracy in the spin temperature/covering fac- 
tor ratio (Tspin//), for the DLAs where the column 
density is known. Therefore, since there is a mix of 
detections and non-detections in DLAs at Zabs ~ 2, 
whereas at Zabs ~ 2 H i 21-cm tends not to be detected, 
iKanekar fc Chengaluil ([2003') suggest an evolution in the 
spin temperature of the gas, where this is exclusively 
high at high redshift. This has direct bearing on the rel- 

* Optical spectra from QSOs come from continuum sources ^ 1 
pc in diameter, whereas radio sources can exceed ~ 100 pc (see 
Imifc et al. 2005). 

^ Which applies t o all but one of the known 21-cm absorbers 
dRoberts et al.|[T97^) . 
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ative proportions of the cold and warm neutral media 
(e.g. iKanekar k. Chengaluil [2003t ISrianand et al.ll2012l ) 
and thus the reservoir of raw material available for star 
formation in the early Universe. 

However, the conclusion that the spin temperature 
evolves with redshift relies on the assumption of a cover- 
ing factor for each DLA, and ICurran et al.l ()2005D have 
shown that, if the degeneracy is left intact, Tspin/f is not 
the order of magnitude larg er at Zabs ^ 2, as suggeste d 
by the Tgpin d istribution of Kanckar fc Chengaluii ( 2003[ ) . 
Furthermore, ICurran fc Webb ( 2006) suggest that geom- 
etry effects introduced by an expanding Universe could 
contribute to the distribution, indicating that the cov- 
ering factor may actually be the dominant parameter. 
Since then there have been further high re dshift searches 
(jCurran et al.l 120101: ISrianand et"all 120121) ^ . This takes 
the number of confirmed DLAs and sub-DLAs searched 
for 21-cm to 25 at z > 2, the same number as at z < 2, 
thus allowing us to adequately address the spin temper- 
ature versus covering factor issue. 

2. NEW Hi 21-CM ABSORPTION SEARCH RESULTS 

Although there are now equal numbers of 21-cm 
searches above and below Zabs — 2, the detection rate 
at Zabs > 2 remains considerably lower (5 out of 25) than 
for the lower redshift sample (14 out of 25). Since the ve- 
locity integrated optical depth of the 21-cm absorption, 
normalised by the total neutral hydrogen column den- 
sity, is proportional to //Tspin (Sect. [T|), this suggests 
that either Tgpin is generally higher at Zabs > 2 or that 
/ is generally lower, or a combination of both, elevating 
the ratio at high redshift. 

In order to obtain limits to the line strengths for the 
non-detections we must assign a profile width to each ab- 
sorber. This is usually assumed, but here we estimate the 
full- width half maximum (FWHM) of the profile from the 
21-cm FWHM— Mg II 2796 A equivalent width correla- 
tion for the detections, specifically FWHM [km s""'^] w 
13 X Wj^2796 (ICurran et al.l l2007bD . Where the 
Mgii equivalent widths are not available (generally at 
Zabs ~ 2.2 where the transition is redshifted out of 
the optical band), we estimate the putative FWHM 
via the metallicity — Mg il equivalent width correlation 
(jMuroh v ct al. 20 07|), sp ecifically W.^^rge ^ 2.0 [M/H] + 
4.0 (jCurran et al.ll2007b() . When neither the Mg 11 equiv- 
alent width nor metallicity are available, we assign a 
FWHM of 25 km s"\ the mean of the detections. 

Given the FWHM of the putative line, in order to ob- 
tain the best signal-to-noise ratio, based upon the de- 
tection of a single channel we would smooth the data 
to a spectral resolution of Av =FWHM. This has the 
effect of improving/degrading the measured r.m.s noise, 
and thus the optical depth, by a factor of Au/FWHM. 
The newly derived integrated optical depth of the line is 
thus 
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° From the Sloan Digital Sky Survey DLA Data Releases 1-3 
|Prochaska ct al. 2005) and the SDSS DR7 DLA Data Release 7 
lINoterdaeme ct al.,,2009.) . respectively. 

Excluding the FWHM = 240 km s"^ outlier d ue to the ab- 
sorber in the complex sight-line towards 3C 336 ICurran et al.l 
I2007al) . 
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Fig. 1. — The lower limits to the spin temperature/covering 
factor degeneracy obtained from the recalculated velocity inte- 
grated optical depths versus those obtained from the previously 
published values. The circles show the values obtained by estimat- 
ing the FWHM from the Mgll equivalent width, the squares using 
a FWHM of 10 km s~^ and the triangles using the mean 25 km s~ 
of the detections. This is the default FWHM used when the Mg II 
equivalent width or metallicity is not available, thus giving the 
overlap between circles and triangles. The dotted line shows where 
the new and previous values are equal. 



= Tpeak X VAwFWHM, (2) 

where Tpoak is the peak optical depth of the putative 
line. This has the overall effect of scaling the previously 
published integrated optical depths, where these are cal- 
culated per channel, by \ f rdv] I \ f rdv] 

^ ' ^ LJ J now ' LJ J previous 



Tpoak X Av 



Tpoak X VAv FWHM _ /FWHM 



Av 



(3) 



We therefore recalculate the limits from the pub- 
lished r.m.s. noise levels and spectral resolutions and 
in Fig. [1] show the effect this has on the Tspin/f lim- 
its compared to the previous values. For the method 
prescribed above (shown as circles), the largest change 
is for the Zabs = 0.5579 absorber towards PKS0118- 
272, for which the limit chm bs from Tspin/f > 850 K 
(|Kanekar fc Chengaluii [200 ID to > 1900 K and, on av- 
erage, the scaled values are 160 K "warmer" than those 
published. As seen from the figure, however, the new 
hmits are mostly concentrated close the [/ Tdv]^^^ = 
\ f rdv] . line, with equal numbers of points to ei- 

J previous 

t her side. 

ICurran et al.l (|2007b[ ) attribute the correlation between 
the velocity spreads of the 21-cm and Mg 11 profiles to the 
neutral and singly ionised gas being spatially coincident 
in 21-cm detected DLAs. However, given the possibility 
that much of the ionised magnesium could arise in warm 
photoionised gas, and is thus not tracing the same sight- 
line as the 21-cm, in Fig. [T] we also show the Tspin/f 
limits derived from assuming a single FWHM. For 10 
km s~^ (e.g. ISrianand et al.|[2012|) . we find the largest 
difference to be 4300 K, with a mean rise of 600 K and 
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for 25 km s^^, the largest difference is 2500 K, with a 
mean decrease of 120 K. Thus, scahngs implementing the 
mean FWHM of the detections have less of an effect than 
using 10 km s~^, although either method is preferable to 
the inhomogeneously calculated published values, since 
each is consistent in smoothing the data to the FWHM 
and thus calculating the limit per smoothed channel. 

In Fig. [2] we show the newly derived limits against 
the absorption redshift. Applying a survival analy- 
sis to th e 3(7 lower limits in Tgpin/ f, via the ASURV 
package (jlsobe et al.l |1986[ ) , we obtain mean values of 
logio(r«pin//), < 2 = 3.25±0.15 and logio(T,pi„//)^ > ^ = 
3.56±0.17. That is, the mean value of Tgpin// at Zabs > 2 
is double that at Zabs < 2 (3600 cf. 1800 K). This 
is considerably less than the ratio of 8.0 between the 
two redshift regimes (8100 cf. 1020 K) found by ap- 
plying the survival an a lysis to the spin temperatures of 
iKanekar fc Ghengaluri ^OOf). Since our recalculation of 
the limits has actually slightly increased the mean value 
of the spin temperature/covering factor degeneracy, we 
attribute this significant difference to the covering factor 
being left as a free parameter. 

3. SPIN TEMPERATURE OR COVERING FACTOR 

The addition of the new data, together with our treat- 
ment of the non-detections, indicates that a factor of 
only « 2 in the spin temperature/covering factor ratio 
between the two redshift regimes need be accounted for. 
However, it remains unclear whether Tspin or / is the 
dominant factor. 

There have been several attempts to quantify the cov- 
ering factor through Very Long Baseline Array (VLBA) 
imaging of the background quasar, in which / is esti- 
mated as the ratio of the compact unr esolved compo- 
nent' s flux to the total radio flux (e.g. iBriggs fc Wolfd 
119831 IKanekar et al.lj2009t ISrianand et al.ll2012D . From 
this. IKanekar et al.l ( 2009 ) conclude that for the 21-cm 
searches, / > 0.4 and hence low covering factors cannot 
explain the dearth of detections at high redshift. How- 
ever, this method gives no information on the extent of 
the absorber nor on how effectively it covers the emission. 
Neither does it give information on the depth of the line 
when the extended continuum emission is resolved out. 

So although an evolution in spin temperature with red- 
shift is implied, this method invokes several crucial as- 
sumptions which could even lead to unphysical results 
(e.g. the - pc obtained for the ex tent of PKS 0405-331 
at 398 MHz, IKanekar et"alll2009| ). On the other hand, 
due to the effects of an expanding Universe, objects at 
1 ^ z ^ A are always close to a single angular diame- 
ter distance (~ 1500 Mpc), irrespective of redshift (e.g. 
iHoggl [T999| ) . Beyond z ~ 1, the angular diameter dis- 
tance decreases with increasing redshift, which will lead 
to a natural imbalance between the low and high-z DLA 
samples. In the small angle approximation the covering 
factor can be defined as 
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dabs -DAqso \ ' 
C^QSO-D^DLA/ 



(4) 



where Oqso is the angle subtended by the background 
emission and o?abs a-nd c?qso the linear sizes of the cool 
absorbing gas and background 21-cm emission, respec- 
tively (Fig. [g. Since absorbers at z > 1 are always at 



a similar angular diameter distance as the background 
quasar (DAdla ~ DAqso), these will therefore, for a 
given system, always have a lowe r coverage than a coun- 
terpart with DAdla < -D^Qso (ICurran fc Webbl l20061. 
an inequality which can only exist at z < 1. 

In Fig. |4]we show the angular diameter distance ratio, 
-DAdlaZ-D^qso, versus the absorption redshift for the 
21-cm searches. Note that this distribution is free from 
any assumptions (including our own with regard to the 
FWHMs of the non-detections). From this we see a mix 
of IJAdlaZ-D^qso at low redshift (the low values are due 
to z < 1 DLAs towards much higher redshift quasars), 
whereas values at z 1 are exclusively high. That is, 
a reproduction of the " spin t emperature" distribution of 
IKanekar fc Chengaluil ()2003l). At least the pattern if not 
the magnitude, given that IKanekar fc Chengalud (|2003l ) 
suggest that the average spin temperature at Zabs > 2 is 
eight times that at Zabs < 2, which is clearly not repro- 
duced in Fig. m 

Referring to Sect. however, we only require TLpin// 
at Zabs > 2 to be about double that at Zabs < 2 in order to 
account for the distribution. From Equ.jH for a given ab- 
sorber and emitter size, / oc 1/(13 Adla/£'^qso)^7 and 
we find mean values of (£'Adla/-D^qso)^2 < 2 — 0.68 

and (-DAdla/-D^qso)^2>2 ~ ^-^(^i a ratio of 1.6. It is 
therefore apparent that the distribution in angular di- 
ameter distance ratios, and their effect on the covering 
factor, is close to that required to account for the ob- 
served Tspin// distribution. 

Updating the statistics with the new results, we find a 
binomial probability of 5.52 x 10~^ of obtaining at least 
9 out of 11 detections at UAabs/^'^QSO < 0-8 and 29 
out of 39 non-detections at ZJ^absZ-D^qso > 0.8, by 
chance.^ That is, a significance of 4.03cr, lending weight 
to the argument that large angular diameter distance 
ratios are correlated with the non-detection of 21-cm ab- 
sorption. Since there is no physical link between the spin 
temperature and angular diameter distance ratio, this 
suggests that, in the absence of any knowledge of the ab- 
sorption cross-section in relation to the emitter extent, 
the covering factor plays a much more major role than 
the spin temperature in the detection of 21-cm absorp- 
tion. 

4. DISCUSSION AND SUMMARY 

With the advent of 21-cm absorption searches in high 
redshift DLAs there has been much debate on whether an 
evolution in spin temperature or lower covering factors 
are responsible for the relatively low number of detec- 
tions at Zabs > 2. New searches of SDSS DLAs, which oc- 
cult the sight-lines to radio-loud quasars, take the num- 
ber searched at these redshifts to the same as at Zabs < 2, 
giving much needed data with which to address this is- 
sue. By giving the 21-cm non-detections a more thorough 
analysis than in other studies, through estimating the 
FWHMs of the non-detected profiles from the observed 
optical properties and applying a survival analysis to the 



^ D A^'as/D^OSO = 0.8 is the cut used by ICurran fc WebbI 
l|2006l ) and, although arbitrary, this is the lowest ratio which gives 
an appreciable enough sample size in the lower bin. When a ratio 
of 0.9 is used there are 38 DLAs with DA^^^^^J DAq^q > 0.9 and 
12 with DA^\,s/DAqqo < 0.9, which increases the significance of 
the correlation to 4.31(t. 
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Fig. 2. — The spin temperature/covering factor degeneracy versus the absorption redshift for DLAs searched in 21-cm absorption. 
The filled symbols/histogram represent the detections, the unfilled symbols/histogram the non-detections (3o" limits), with the shapes 
representing the type of galaxy with which the DLA is associated: circle-unknown type, star- spiral , square-dwarf, triangle-low surface 
brightness (LSB). The coloured symbols show the recent high redshift searches of SDSS DLAs llCurran et al.ll2010l : ISrianand et al.ll20l2l 'l. 
In the histogram the limits to T^pi^/ f are treated as actual values. 




Fig. 3. — Absorber and quasar cross-sections with respect to the 
angular diameter distances. If the angle required to fully subtend 
d^ljs is Iess than that to subtend dqso then / < 1. 



lower limits in Tspm/ f, we find: 

1. That mean the spin temperature/covering factor 
ratio at Zabs > 2 is only double that found at 
Zabs < 2. This is significantly lower than the fac- 
tor of eight in spin temperature, determined from 
previous estimates of the covering factor. 

2. This reduced factor can be accounted for by the 
effects of an expanding Universe, where the angular 
diameter distance at Zabs > 1 is always close to, and 
sometimes exceeds, that of the background radio 
source. That is, the high redshift DLAs are always 
disadvantaged in how effectively they can intercept 
the background radio flux for a given DLA-QSO 
pair. 

3. Since the non-detections tend to be associated 
with high ratios in the DLA-QSO angular di- 
ameter distances {DAahs/DAQso ~ 1), whereas 
the detections are associated with low ratios 



{D AgCbs / D Aqso < 1), at aZ^ redshifts, this strongly 
suggests that the coverage of the background radio 
source by the absorber is the dominant factor. 

This runs contrary to previous studies which estimate 
the covering factor via VLBA imaging, which support 
an evolution in the spin temperature. However, these 
estimates are based upon several pivotal assumptions, 
whereas the angular diameter distance ratios are not. 

We therefore suggest that there is no evidence for 
the spin temperature to increase with redshift, which 
implies no evolution in th e fraction of t h e warm neutral 
gas content of DLAs. ICurran et al.l ()2007b( ) find a 
rspin//~iiietallicity anti-correlation. If the spin tem- 
perature does not evolve with redshift , as does th e 
metallicity (Prochaska et al. 2003; Curr an et al.l 120041 ) , 
this may imply that the metallicity is proportional to 
the covering factor, which could be due to the larger 
absorbing cross-sections of the larger galaxies, as evident 
through the spirals which have t he largest me tallic- 
ities and the lower values of Tgpin/ f (jCurran et alU lOlO). 

This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and 
Space Administration. This research has also made use 
of NASA's Astrophysics Data System Bibliog raphic Ser- 
vice and ASURV Rev 1.2 (jLavallev et al.lll992h. which irn - 
plements the methods presented in llsobe et al.l (|1986l ) . 
The Centre for All-sky Astrophysics is an Australian 
Research Council Centre of Excellence, funded by grant 
CE11E0090. 
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